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We have analyzed the latest NA50 data on J/tp production 
in pA and AA collisions. The J/tp production is assumed 
to be a two step process, (i) formation of cc pairs, pertur- 
batively calculable, and (ii) formation of J/tp from the pair, 
a non-perturbative process, which is conveniently parameter- 
ized. In a nuclear medium, as the cc pair passes through the 
nuclear medium, it gain relative square momentum and some 
of the pairs can gain enough square momentum to cross the 
threshold for open charm meson, leading to suppression in nu- 
clear medium. Few parameters of the model were fixed from 
the latest high statistics NA50 pA and NA38 SU total J/ip 
cross sectional data. The model then reproduces the central- 
ity dependence of J/ip over Drell-Yan ration in 200 GeV/c 
S+U and 158 GeV/c Pb-hPb collisions. 



I. INTRODUCTION 

Recently, in Quark Matter 2002, NA50 collaboration 
presented their analysis of the nuclear absorption of J/ip 
in high statistics 450 GeV pA collisions [1]. They esti- 
mated the J/ip nucleon absorption cross section {o'^^^) 
in the framework of Glauber model. High statistics 450 

4.4 ± 1.0 mb [1]. They also 



GeV pA data yield aH'^^ 



estimate a common cr^/^^ from latest pA and NA38 200 

GeV/c S-KU data [2], a^l^^ =4.4 ± 0.5 mb. The ex- 
tracted absorption cross section is much smaller than the 
earlier value of 6.4 ± 0.8 mb extracted from fit to ear- 
lier NA5G data [3] or 7.1 ± 3.0 mb obtained from a fit to 
NA38 S-f U data [2]. Within error, the S-l-U cross sections 
are compatible with pA cross sections. In Quark matter 
2002, NA50 collaboration also presented their prelimi- 
nary analysis of 2000 Pb+Pb run [4] . Et dependence of 
the J/ip over Drell-Yan ratio, in the rapidity range of 2.9- 
4.5 was presented. Compared to 1998 run [5], 2000 data 
are flatter, suppression being more at low Et and less at 
high Et- The difference may be attributed to the differ- 
ence in method of analysis. In 1998, NA50 collaboration 
presented combined result of standard as well minimum 
bias analysis. In the 2000 run, only the standard analysis 
was performed. Centrality dependence of J/t/i absorption 
is still anomalous in the sense that Glauber model of nu- 
clear absorption fails to explain the data. The 1998 data 
of NA50 collaboration [5] were analyzed in a variety of 
models, with and without assumption of QGP [6-9,11]. 
We have shown that a QCD based model, where J/-0's 



are absorbed in nuclear medium could explain the data 
[9,11,12]. Parameters of the model were obtain from a fit 
to then existing pp/pA/AA total J/ip cross section data 
[3]. 

Preliminary analysis of high statistics pA data [1] 
with the implication of less absorption of J /ip in nu- 
clear medium, along with the changed suppression pat- 
tern observed in the preliminary analysis of NA50 2000 
Pb-fPb run [4], necessitated the re-examination of the 
QCD based nuclear absorption model [9]. In the present 
brief report, we have analyzed the latest NA50 data on 
pA/SU total J I'll) cross section [1] to obtain the param- 
eters of the QCD based nuclear absorption model. The 
model was then used to explain the centrality dependence 
of J /ip over Drell-Yan ratio in 200 GeV/c S+U collisions 
and in 158 GeV/c Pb-|-Pb collisions (preliminary data). 
As will be shown in the following, once the model param- 
eters are fixed from J /ip total cross sectional data, the 
QCD based nuclear absorption model give a consistent 
description of the centrality dependence of J /ip suppres- 
sion in S-l-U collisions as well as in Pb+Pb collisions. The 
plan of the paper is as follows: in section 2, we will briefly 
describe the model and obtained the model parameters 
from the new high statistics pA/AA NA50 data. In sec- 
tion 3, the NA50 data on Et dependence of J /ip over 
Drell-Yan ratio in S+U and in Pb+Pb collisions are an- 
alyzed. Summary and conclusions are drawn in section 
4. 



II. J/iP PRODUCTION AND SUPPRESSION 

In the QCD based nuclear absorption model [9,12], 
J/'ip production is assumed to be a two step process, 
(a) formation of a cc pair, which is accurately calculable 
in QCD and (b) formation of a J /ip meson from the cc 
pair, which is non-perturbative but can be conveniently 
parameterized. The J /ip cross section in AB collisions, 
at center of mass energy -y/s is written as, 
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where runs over all parton flavors, and — + 

4rn^. The K factor takes into account the higher order 
corrections. The incoming parton momentum fractions 
are fixed by kinematics and are Xa — (V ^'f + 4Q^/ s + 
x_f)/2 and Xb = i\/ Xp + 4Q^/s - xf)/2. aab-*cc are 
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the sub process cross section and arc given in [14]. 
Fcc^j/^{q^) is the transition probabihty that a cc pair 
with relative momentum square evolve into a physical 
J/tjj meson. It is parameterized as, 
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FIG. 1. The experimental ratio of total J/tp cross section 
and Drell-Yan cross sections in proton-proton, proton- nucleus 

and nucleus-nucleus collisions. The fit to the data obtained 
in the QCD based nuclear absorption model is shown as solid 
lines. 

In a nucleon- nucleus/nucleus- nucleus collision, the pro- 
duced cc pairs interact with nuclear medium before they 
exit. It is argued [12] that the interaction of a cc pair with 
nuclear environment increases the square of the relative 
momentum between the cc pair. As a result, some of 
the cc pairs can gain enough relative square momentum 
to cross the threshold to become an open charm meson. 
Consequently, the cross section for J /ij.) production is re- 
duced in comparison with nucleon-nucleon cross section. 
If the J meson travel a distance L, q^ in the transition 



probability is replaced to q 



e L, e being the 



relative square momentum gain per unit length. Param- 
eters of the model {ap^KNj/^ and e^) can be fixed from 
experimental data on total J /ip cross section in pA/AA 
collisions. Earlier NA50 data [3] are well described with 
KNj/^^QAb%, = 0.225GeV'^/fm and ap = 1.0 [9]. 
As mentioned in the beginning, NA50 collaboration pre- 
sented high statistics pA data on J/tp cross section. They 
have measured B f^^^cr (J / ip) / (7{DY) . In Fig.l, experimen- 
tal data are shown as a function of nuclear length. The 
Drell-Yan cross sections donot have any A or alternately 
any L-dcpcndence. The observed L dependence is then 
due to J ftp's only. We fit the data with two parameters, 
Nnorm = KNj/pgi/a^^{nb) and square momentum gain 
factor (ap being kept fixed at 1). In Fig.l, the fit ob- 
tained to the data is shown. The two sets of data at 200 
GeV/c and 450 GeV/c could be fitted with a common 



square momentum gain factor, e^=: 0.1875 GeV"^ / fm, a 
value 20% lower than the value obtained earlier [9] . Low- 
ering of indicate less absorption of J/i/j'a in nuclear 
medium, in agreement with the Glauber model calcula- 
tions. While the square momentum gain factor do not 
show energy dependence, the evident energy dependence 
of the cross section ratios shows up in the other param- 
eter of the model Nnorm- We obtain Nnorm = 10.18 at 
200 GcV/c and iV„o„„ = 4.43 at 450 GcV/c. The energy 
dependence of J ftp cross section being taken care of in 
the model (Eq.l), the energy dependence of Nnorm is due 
to the Drell-Yan cross sections only. In the mass range, 
2.9 > M > 4.5 GeV, the Craigie parameterization [15], of 
the Drell-Yan cross section, a{DY) oc e~^^-^^/ gives 
for the ratio (T(-Dy)450Gey/o'(-DF)2ooGeV' = 2.1 -3.1, con- 
sistent with the presently obtained ratio of 2.29. 



III. Et dependence of j/^p/dy ratio 

In the present section, we analyzed the Et dependence 
of J ftp over Drell-Yan ratio in 200 GeV/c S-l-U collisions 
and in 158 GeV/c Pb-I-Pb collisions. Details of the cal- 
culation can be found in ref. [9]. 
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FIG. 2. The transverse energy dependence of J/tp over 
Drell-Yan ratio in 200 GeV/c S+U collisions. The dashed 
and dotted lines are the fit obtained in the Glauber model 
of nuclear absorption with cr^^*^=7.1mb and 4.4 mb respec- 
tively. The solid line is the fit obtained in the QCD based 
nuclear absorption model. 



Fig. 2, NA38 experimental ra- 

^^a{J/ip)/a{DY) for 200 GeV/c S-fU collisions 
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[2] is compared with the present model calculation (solid 
line). The data are well described in the model. For 
comparison purpose, we have also shown the fit obtained 

to the data in the Glauber model. The dashed line is the 



Glauber model calculations with a 



J/tljN_ 



7.1 mb. How- 



ever, quality of fit gets poorer if we use lower value (4.4 

mb) for (y'^lf'^ as cixtractcd from the most recent data 
(the dotted line). It is evident that in the Glauber model 



2 



of nuclear absorption, the transverse energy dependence 
of J/ij/'a prefer larger absorption of J/tp in S+U collisions 
than indicated in pA/AA total J/ip cross section data. 
The QCD based nuclear absorption model, on the other 
hand, give consistent description of the data. 
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FIG. 3. NA50 preliminary transverse energy dependence 

of J/tp over Drell-Yan ratio in 158 GoV/c Pb+Pb collisions. 
The dotted line is the Glauber model calculation. The dashed 
line is the QCD based nuclear absorption model calculation 
without including the effect of Et fluctuations. The solid line 
is obtained with including the effect of Et fluctuations. 

Transverse energy dependence 

of the ratio, B ^^a{J / 'tp) / a{DY) in 158 GeV/c Pb+Pb 
collisions could not be explained in the Glauber model 
of nuclear absorption. The preliminary NA50 data [4]ob- 
tained in 2000, are flatter compared to 1996/1998 data 
[5], still the usual Glauber model of nuclear absorption 
fails to describe the data. In Fig. 3 the experimental data 
are shown along with the Glauber model fit to it obtained 
with (T^/^ =4.4 mb. The Glauber model calculation (the 
dotted line) agree only for the peripheral collisions. For 
more central collisions, the model predict less absorp- 
tion than in experiment. Also the Glauber model shows 
saturation at more central collisions, in contrast to ex- 
periment. The 2nd drop in the ratio is mainly due to 
transverse energy fluctuations at a fixed impact param- 
eter [6,10]. In the present model, we take into account 
Et fluctuations at a fixed impact parameter h, by the 
replacement: L{h, s) L{b, s)Et/ < Et > (b). In Fig.3, 
the dashed is obtained in the QCD based nuclear absorp- 
tion model, without incorporating the Et fluctuation ef- 
fect. It explains the data up to 100 GeV, the knee of 
the Et distribution, beyond which only, Et fluctuations 
are important. U Et fluctuation is included, the data are 
explained throughout the Et range (the solid line). Con- 
sistent description of Et dependence of J/ tp over Drell- 
Yan ratio in 200 GeV/c S-hU and 158 GeV/c Pb+Pb 
collisions clearly demonstrate that the QCD based nu- 
clear absorption model, together with transverse energy 
fluctuation could explain consistently the J/ %p absorption 
data at SPS energy. 



IV. SUMMARY AND CONCLUSIONS 

To summarize, we have analyzed the latest NA50 data 
[4] on J/ip suppression obtained in 2000 Pb+Pb run. The 
data were analyzed in the QCD based nuclear absorption 
model. The J/V' production is assumed to be two step 
process, (i) production of cc pair, purturbatively calcu- 
lable and (ii) formation of J/tp from the cc pair. In a 
nuclear medium, due to random collisions with medium, 
relative square momentum of the cc increases and some 
of the pair may gain enough square momentum to cross 
the threshold for open charm meson. We obtain the pa- 
rameter, square momentum gain factor per unit length, 
e^, from a fit to the latest NA50 data [1], and obtain 
£^=0.1875 GeV"^ / fm. The value is lower than the value 
e^= 0.225 GeV"^/ fm, obtained in our earlier work [9]. 
There, we had fitted the earlier version of the NA50 data 
[3]. Reduction in indicate that the J/ip^s are less ab- 
sorbed in nuclear medium. This is in accordance with 
the J/'^-nucleon absorption cross section, extracted from 

the data [1,3]. o'l^lf^ extracted from the earlier version 

of the NA50 data [3] was large, a^^*^ = 6.4 ± 0.8 mb, 
compared to the value extracted from the latest data set 
[1], o-'^l^^ = 4.4 ± 1.0 mb. With the new model param- 
eter, the model could explain the centrality dependence 
of J ftp over Drell-Yan ratio in 200 GeV/c S+U collisions 
as well as the latest NA50 data on the centrality depen- 
dence of J/tp over Drell-Yan ratio in 158 GeV/c Pb+Pb 
collisions. 



* e-mail:akc@veccal.ernet.in 
[1] NA50 collaboration, P. Cortese et al, in Quark Matter 

2002, Nantes, Prance (unpublished). 
[2] NA38 Collaboration, M. C. Abreu at. al., Phys. Lett. 

B449, 128 (1999). 
[3] NA50 collaboration, M.C. Abreu et. al., Phys. Lett. 

B410, 337 (1997). 
[4] NA50 collaboration, L. Ramello et al, in Quark Matter 

2002, Nantes, Prance (unpublished). 
[5] NA50 collaboration, M. C. Abreu et al. Phys. Lett. B 

477,28(2000). 

[6] J. P. Blaizot, P. M. Dinh and J.Y. Ollitrault, Phys. Rev. 

Lett. 85,4012(2000). 
[7] A. K. Chaudhuri, Phys.Rev. C64,054903(2001), Phys. 

Lett. B527,80(2002) 
[8] A. Capella, E. G. Ferreiro and A. B. Kaidalov, hep- 

ph/0002300, Phys. Rev. Lett. 85,2080 (2000). 
[9] A. K. Chaudhuri, Phys. Rev. Lett. 88, 232302(2002). 
[10] A. K. Chaudhuri, Phys. Rcv.C66,021902 (2002). 
[11] A. K. Chaudhuri, nucl-th/0207082. 

[12] J. Qiu, J. P. Vary and X. Zhang, hep-ph/9809442, Nucl. 
Phys. A698, 57wl (2002). 



3 



[13] NA50 collaboration, M.C. Abreu et al, Phys. Lett. 

B410,337 (1997). 
[14] C. J. Benesh, J. Qiu and J. P. Vary, Phys. Rev. C50, 1015 

(1994). 

[15] N.S. Craigie, Phys. Reports, 47 (1978) 1. 



4 



